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SUMMARY 

An experimental investigation was conducted in a 3.84- by 10- inch 
trannel to determine the mass transfer by siiblimation, heat transfer, 
and skin friction for an iced surface on a flat plate for Mach nvanbers 
of 0.4, 0.6, and 0.8 and pressure altitudes to 30,000 feet. Measure- 
ments of rates of sublimation were also made for a Mach number of 1.3 
at a pressure altitude of 30,000 feet. 

The results show that the parameters of sublimation and heat trans- 
fer were 40 to 50 percent greater for an iced surface than was the bare- 
plate heat -transfer parameter. For iced surfaces of equivalent rough- 
ness, the ratio of sublimation to heat-transfer parameters was found to 
be 0.90. The sublimation data obtained at a Mach n\miber of 1.3 showed 
no appreciable deviation from that obtained at subsonic speeds. The 
data obtained indicate that sublimation as a means of removing ice 
formations of appreciable thickness is usually too slow to be oi much 
value in the de-icing of aircraft at high altitudes. 


INTRODUCTION 

Study in the field of icing meteorology has indicated that the 
greatest probability and frequency of aircraft icing encoimters exist 
at altitudes below 20,000 feet (ref. l). Modern high-speed, high- 
altitude aircraft might thus be expected to climb rapidly through icing 
conditions with subsequent removal of residual ice formations outside 
the icing conditions. In recent years there has been speculation con- 
cerning the possibility of removing these residml ice formations during 
flight in clear air at hi^ altitudes by the process of sublimation 
(mass transfer from solid to vapor), thereby minimizing the icing- 
protection requirements for aircraft. Similarity between the processes 
of sublimation of ice and evaporation of water makes the principles 
evolved by previous investigators in the field of evaporation applicable 
to sublimation studies. From Reynolds analogy and from the work of many 
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investigators, notably ColbTirn (refs. 2 and 3), in correlating experi- 
mental data, the relations between mass transfer, heat transfer, and 
skin friction for a flat plate are well known. Thus a reasonably accurate 
prediction of the rate of sublimation of a smooth ice surface subjected 
to a tangential low-velocity air stream could be made from available 
data. However, natural ice formations resulting from the impingement 
and freezing of small water droplets on a surface are rarely smooth, and 
very little is known of the sublimation rate of roijgh ice in a high- 
velocity air stream. 

The present investigation was \indertaken to determine the rate of 
sublimation from the iced surface of a flat plate in a subsonic and a 
low supersonic air stream, and to relate the siiblimation and heat trans- 
fer from a roiagh ice surface. The investigation was conducted in a 
3.84- by 10-inch t\mnel at the NACA Lewis laboratory. 


SmBOLS 

The following symbols are used in this report: 

Cf average skin-friction coefficient 

Cf local skin-friction coefficient 

Cp specific heat at constant pressure, Btu/(lb)(°F) 

e vapor pressure, in. Hg abs 

g acceleration due to gravity, ft/sec^ 

hg unit thermal convective conductance, Btu/(hr)(sq ft)(*^/ft) 

hjjj unit mass-transfer conductance, lb/(hr)(sq ft)(lb/lb) 

k roughness projection height (ref. 4) 

Lg latent heat of evaporation, Btu/lb 

Lg latent heat of sublimation, Btu/lb 

M Mach number 

m^^ molecular weight of air 


“w 


molecular weight of water 
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Pr^ Prandtl nvunber for heat transfer 

h 

PTjjj Prandtl nvunber for mass transfer 

p static pressure, in. Hg abs 

q rate of heat transfer per unit area, Btu/(sec)(sq ft) 

Re Reynolds number, pUx/gp 

Stv, Stanton number of heat transfer, - 75 — rr^‘ z — r 

^ PUCp(t^d - tg) 

Sts Stanton number of sublimation, 

t temperature, ^ 

tad adiabatic wall temperatxore, ^ 

U velocity of fr^ stream, ft/sec 

u velocity component in boundary layer parallel to surface, ft/sec 
W , sublimation rate, lb/(sec)(sq ft) 

X distance from leading edge of model, ft 
y coordinate normal to siirface, ft 

X ratio of specific heats 

5 boundary-layer thickness, ft 

6 * boundary- layer displacement thickness, ft 

6 boundary- layer momentum thickness, ft 

\i viscosity of air, (lb)(sec)/sq ft 

p weight density of air, Ib/cu ft 

Subscripts : 
s surface 

y at distance y above surface 
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0 free stream 

1 outer edge of boundary layer 


ANALYSIS 

The determination of the rate at which ice will sublimate from a 
s\irface in a high-velocity air stream is dependent primarily on the 
knowledge of the absolute value of the dimensionless parameter of mass 
transfer, herein referred to as the Stanton number of sublimation. In 
order to determine experimental values of the Stanton number of subli- 
mation, the method used in reference 5 for the evaporation of water has 
been applied to the sublimation process. By this method the tenipera- 
ture of a wet surface in an air stream ce(,n be calculated; and if the 
quantity of ma.tter evaporated can be meas\ired, experimental values of 
the dimensionless parameter of mass transfer can be determined, pro- 
vided the relation between the heat transfer and mass transfer from the 
surface is known. This relation is cited in reference 6 as the ColbTirn 
relation and is given as (in the terminology of this report) 



hmCp \Pr^/ 


The Prandtl nvimber of mass transfer Prj^ is perhaps better known as 

the Schmidt number or Taylor number. The Lewis relation also cited ’in 
reference 6 is a special case of the Colburn relation in which thd ratio 
of the Prandtl numbers is equal to unity. For the case of water evapo- 
rating into air, the value of the preceding relation is, nearly eqml to 
unity (refs.' 5 and 7) and the ratio of the dimensionless parameters 
(Stanton numbers) of heat and mass transfer is also nearly eqxial to 
unity. 


In the treatment of the problem given herein, the equation for the 
kinetic temperature of a wet surface .from reference Sis used as the 
starting point. In the notation used herein, this eqtiation has the 
form 

t -t f,, 

ad s stj, Cp ' 

Analogy between the processes of evaporation and sublimation makes, 
possible the substitution of the latent heat of sublimation Ls for the 
latent heat of evaporation Lg and the substitution of the Stanton n-um- 
ber of sublimation Stg for the Stanton number of evaporation Stg in 
equation (l), which becomes 
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■'ad 


t = !^ 

s Stj^. maVPl - Pi - ® 


1/ '"p 


( 2 ) 


This eqmtion resiolts from combining the equation for the Stanton number 
of sublimation as defined by 


Sts = 


W 


pu 




\Pi 


(3) 


- e. 


Pi 


“ 


with the equation for the Stanton nuinher of heat transfer 


St^ = 
h 


PUCp(taa - tg) 


(4) 


where 


q = LgW (5) 

In the determination of values for Stg from equation (3), the weight 

rate of sublimation W is determined experimentally and all the other 
terms are either known or directly measxired except the vapor pressure 
at the iced surface 63 . The values of es correspond to the vapor 
pressure over ice at the surface temperature ts. Both ts and 63 
may be obtained from equation ( 2 ) by a trial-and-error process provided 
the ratio Stg/Stj^ is known. The Stanton number ratio has been re- 
lated to the two-thirds power of the ratio of Schmidt n\miber to Prandtl 
nioiriber;, and as both of these numbers are independent of the initial state 
(solid or liquid) of the material transferred, there appears to be no 
reason to differentiate between sublimation from an ice siirface and 
evaporation from a wetted surface. 

By means of eq'uations (2), (3), and (i) it is possible to compute 
values of Stj^ and Stg based on measured values of W and q. Be- 
cause a true measxxre of tg, the temperature of the ice-air interface 

at which sublimation is taking place, is difficult to obtain, it is 
necessary to compute its value along with the value of eg, the equi- 
libritun vapor press\ire over ice corresponding to tg. The psychrometric 
tables relating tg and eg and the three equations (2), (3), and (4). 

permit the evalmtion of the four unknown quantities tg, eg, Stj^, 
and Stg. 
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Algebraic diff icxilties require that an iterative process be em- 
ployed for the evaluation of the four unknowns. This iterative process 
begins with the use of equation (2) along with the psychrometric tables 
to compute a first approximation for tg and eg for a given set of 
known experimental values of t^d.^ Pl> e^. In order to make this 

coinputation, it is necessary to use a trial value of Stg/Stj^. As a 
starting point, the value of unity is employed. With tg and es com- 
puted in this way, values of Stg are obtained from equation (3) by 
using the measured values of W corresponding to the experimental 
values of tg^^, p^, and e^^. Experimental values of Stj^ which do not 

depend upon the solution of equation (2) were determined from equation 
( 4 ) by conducting an investigation which was independent of the subli- 
mation study in a manner such that values of the heat flow q, and the 
surface temperatvire tg, could be measured. These computed val\aes of 

Stg and Stjj are then used to form a new ratio Stg/Stj^ to be used 
in equation (2) to obtain a second approximation for tg and es. New 
values of Stg are then computed by means of equation (3). This pro- 
cess is repeated until successive values of the ratio of Stg and Stjj 

differ by an amount that is no longer significant in view of the order 
of accuracy of the experimental measurements. 

Iharing the experimental investigation, differences in the ice 
roughness occttrred during the sublimation process. These differences 
were a function of the Mach mmber at which the ice was formed. In 
order to correlate the sublimation Stanton nmnber and the heat-transfer 
Stanton number and to obtain true values of the ratio of Sts/Sth, sub- 
limation and heat-transfer data must be obtained at equivalent values 
of ice roughness. For this purpose and also for the purpose of deter- 
mining the magnitude of the effect of ice on the skin friction and con- 
sequently on the heat-transfer rate, determinations of the skin friction 
both with and without ice on the surface were made. 

In reference 8 the Karman momentum equation for two-dimensional 
flow from which local skin-friction coefficients can be calculated is 
expressed as 


(2 - 4 . 


The momentum thickness 9 and the displacement thickness 8* are de- 
fined by the integrals 
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The average skin-friction coefficient for a surface is the inte- 
grated value of the local skin-friction coefficients; thus 


= 


Pc dx - i P /— + e 

(2 - + ^*!q 


io " I"" 

3(1 

dx / 


dx 


( 9 ) 


The second term inside the integral on the right-hand side of equa- 
tion (9) represents the contribution to the skin-friction coefficient 
due to the pressure gradient along the model. For the negative pres- 
s\ire gradient which exists in a tunnel of constant cross-sectional area, 
the pres sure -gradient term of eqmtion (9) is positive; hence its in- 
clusion in the calculations would have the effect of increasing the 
average skin-friction coefficient, especially at the higher Mach numbers 
(dM/dx increases with increasing test Mach number). For this investi- 
gation the pres sure -gradient term was not included because of the dif- 
ficulty of making the necessary boundary-layer measurements with the 
equipment used for the sublimation work. With the effect of pressure 
gradient neglected, equation (9) becomes 



from which calculations of the average skin-friction coefficient were 
made. 


APPARATUS 

A schematic diagram of the 3.84- by 10- inch tunnel is presented in 
figure 1 and shows the inlet diffuser section with screens, the plenum 
chamber with flow- straightening tubes, the bellmouth tunnel entry, the 
test section, and the outlet diffuser section. The tunnel is designed 
to provide a range of subsonic Mach numbers from 0.3 to 0.8 and a super- 
sonic Mach number of 1.3 over the model. 
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A supply of refrigerated air initially at approximately -20° F and 
with a specific humidity of 5. 0X10““^ pounds of water per pound of dry 
air was conditioned to provide the desired temperatures and humidities 
at the tunnel test section. The humidity of the air stream was con- 
trolled hy means of steam injected. at a point sufficiently far tgjstream 
to ensure thorough mixing at the tunnel entry. 

One wall of the tunnel contained a large glass section for observa- 
tion and visual measurements. The other wall had five portholes for 
access to the inside of the t-unnel and removable plugs for the installa- 
tion of instrumentation at various stations along the tunnel. Permanent 
instrumentation of the tunnel at the test section includes static- 
pressure taps along the top and bottom surfaces of the tunnel and pres- 
suire taps and thermocouples in the plenum chamber. 

The flat-plate model (fig. 2) used for the study was made of wood 
to minimize conduction losses. It was 0.75 inch thick, 3.84 inches 
wide, and 18 inches long. The. formation and sublimation of ice took 
place on a 3.84- by 5.75-inch copper plate set flush in the upper sur- 

7 

face of the model with the leading edge of the copper plate located 4— 

IS 

inches from the leading edge of the model. The upper surface of the 
copper plate was instriamented with fiv'e thermocouples installed flush 

with the surface and spaced at l^-inch intervals along the center line. 

The copper plate ^ hereinafter referred to as the "cold plate," was 
the top surface of a multipass copper box through which cold alcohol 
could be pumped. The alcohol was cooled by pumping it through a coil 
of copper tubing immersed in an alcohol - dry-ice bath contained in a 
Dewar flask. 

Thermocouples located in the inlet and outlet lines to the copper 
box were used in conjunction with a potentiometer of high sensitivity 
to measure the alcohol inlet temperature and alcohol temperature rise 
between the inlet and outlet lines. ' It was found that temperature dis- 
tributions in the inlet and outlet alcohol flow were so great as to 
preclude the use of any sirrple thermocouple installation in the tubes. 

A set of five thermocouples installed in the plane of a cross section 
of the tubes showed that temperature differences of over 20^ F existed 
in the l/4-inch-diameter tubes. The use of Dewar flasks as plenum 
chambers close to the model reduced the temperature distribution range 
in the inlet alcohol flow to less than 0.5^ F. Because temperature 
layers were found to exist in the outlet Dewar flask, a glass stirrer 
was installed in the flask and the temperature of the agitated mixture 
was measured. In order to check the validity of the use of the 
agitated-mixture ten5)erature as the ten^erature at the outlet from the 
copper box, the direction of alcohol flow was reversed and no measurable 
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difference between the temperatures as measured in the flask or in the 
outlet tube was observed. Therefore the tender at ure loss in the line 
between the model and the flask was negligible. A pictorial schematic 
diagram of the alcohol and thermocouple systems is included in fig- 
ure 2. 

Measurement of ice thickness while the tunnel was in operation was 
possible by the use of a short -focal-length telescope mounted on a 
vernier carriage. 

Boundary-layer total -pres sure measurements were made by using a 
16-tube rake mounted as illustrated in figxire 2. The rake tubes were 
0.017 inch in diameter and s/s inch long from tip to rake body. The 
tubes were spaced 0,03 inch between centers for the first l/4 inch 
above the model and 0.04 inch between centers for the remainder of 
the rake. 


PROCEDURE 

The study of sublimation and heat transfer from an iced surface 
requires that the ice surface and the physical characteristics of the 
ice be readily reproducible, and for this investigation it was desired 
that the ice approximate that formed on aircraft surfaces. Preliminary 
experimentation with several methods of forming ice together with lim- 
itations of space and access to the inside of the tunnel indicated that 
ice formed by the condensation of moisture from a humidity-controlled 
air stream on the cold surface of a metal plate most nearly suited the 
requirements . 

The study of the heat-transfer process from the iced surface im- 
poses two additional requirements: (l) the transfer of heat must take 

place without the transfer of mass to avoid interaction of the two 
processes, and (2) the measurement of the temperature of the ice sur- 
face should be possible. The condensation method of forming the ice 
offers an indirect method for obtaining the surface temperature by 
measurement of the humidity of the air stream under conditions for which 
the ice thickness is constant (no net condensation or sublimation 
effect). The ice surface temperature under these conditions is equal 
to the temperature of saturation of the air stream. Thus both the 
additional requirements of the heat-transfer study can be satisfied 
thro"ugh use of the condensation method of forming the ice. Therefore 
for this investigation ice was caused to form on the cold plate by in- 
creasing the humidity of the air stream by means of steam injection. 

When the dew point of the air stream exceeded the surface temperature 
of the cold plate, condensation in the form of a dense frost-like ice 
formation occurred, as shown in figure 3, 
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The movable probe shown in figure 3 was used for purposes of 
focusing "the "telescope and aiding in "the de"termina"tion of ice "thick- 
ness. Throughout the investigation, each ice formation was deposited 
on the refrigerated surface of the cold plate at nearly constant values 
of 5000-foot pressure altitude and 45° F total ten^ierature of the air 
stream in order to assure reproducibility of ice formations for the 
different altitude conditions at each Mach n\miber. During the subli- 
mation study when the desired ice thickness was attained, the ice for- 
mation was caused to sublimate by reducing the humidity of the air 
stream. At the time the humidity of the air was lowered, the desired 
sublimation altitude was set, the alcohol was drained from the copper 
box to minimize the heat-sink effect, and the total temperatiire of the 
air stream was reduced to an arbitrary subfreezing temperature (27° F). 
The sublimation of the ice was in some cases complete in as short a 
period of time as 9 minutes; therefore it was impractical, from the 
standpoint of the time required, to lower .the temperature of the air 
stream to -values of total temperat"ure corresponding to each Mach ntmi- 
ber and altitude condition. The temperature of 27° F was selected be- 
cause it was readily obtained and required approximately the same period 
of time to be established for each of the sublimation conditions. 

Sublimation rates were obtained at Mach numbers of 0.4, 0.6, and 
0.8 for "oressure-altitude conditions of 5000, 15,000, and 30,000 feet. 
Because of test-facility limitations, the Mach number condition of 1.3 
was investigated at a pressure altitude of 30,000 feet only. In order 
to obtain measurements of ice density, a 2- by 4-inch sheet of brass 
0.002 inch thick was wetted on one side and frozen to the copper cold 
plate. Ice was formed on the thin brass sheet by condensation, the ice 
thickness was measiored, and then the brass sheet and the ice were lifted 
from the cold plate and the ice weight and density were determined. 

The ice density was ass\imed to be uniform throughout its thickness and 
at all points on the surface. The rates of sublimation were determined 
from ice-density measurements and plots of ice thickness against time. 
The slope of the secant on the plot of ice thickness against time d"ur- 
ing the interval from 2 minutes after sublimation started to the time 
when sublimation had caused any small portion of the cold plate to be- 
come bare of ice was taken as the rate of change of ice thickness with 
time. The initial 2-minute interval was allowed for stabilization of 
the temperature conditions of the copper box. 

For the heat -transfer study, ice was caused to form on the cold 
plate under the subsonic Mach nuinber and press"ure-altitude conditions 
as previously stated in this section for the sublimation study. After 
the desired thickness of ice had been formed on the cold plate (approx- 
imately 0.2 cm), the humidity of the air stream was reduced to the con- 
dition under which the ice thickness remained constant, indicating that 
the sublimation of the ice and the condensation of the vapor on the ice 
surface were taking place at the same rate, and that the temperat"ure of 
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the ice svirface was the temperatxjre of the frost point of the air 
stream. The frost -point ten 5 )erature was measured from continuous samples 
of the air taken from the tunnel plenum chamber and expanded to cor- 
respond to the static pressxjre existing at the cold plate. A sensitive 
dew-point meter was used for these measurements. In order to provide a 
basis of comparison with previous investigators and for the pxirpose of 
determining the effect that ice on a s\xrface has on the heat -transfer 
rate, data were also obtained for the model with no ice on the surface 
for a similar range of conditions. 

The heat transferred through the ice or to the hare surface was 
determined from the heat gained by the alcohol flowing through the cop- 
per box in the model, as obtained from measurement of alcohol flow rate 
and change in alcohol temperature. The alcohol flow rate was measured 
by diverting the flow by means of the three-way valve shown in figure 2 
and by obtaining volume measvirements of the alcohol over a period of 
45 seconds. The heat gained by the alcohol includes that transferred 
by conduction throijgh the model, which must be subtracted from the total 
to determine the heat gained by convection only. By covering the cold 
plate with a good instilating material of known conductivity, the heat 
gain through this insulating material could be determined and subtracted 
from the total for the insulated model, which leaves the heat gained by 
conduction through the wood. A cork insulator was cemented to the cold- 
plate surface; and the heat gained by conduction through the wood was 
determined for the various Mach numbers, press\ire altitudes, and alcohol 
biilk temperatures. The heat gained by conduction throu^ the cork was 
computed from average plate s\orface temperatiire measurements beneath the 
cork and the adiabatic wall temperature for the stream side of the cork 
together with a value of the conductivity of corkboard from reference 
9. The use of the adiabatic wall temperature as the surface temperature 
of the cork during heat transfer throtigh the cork results in an error 
of the order of 5.0 percent of the heat flow through the cork. However, 
because the heat conducted through the cork is a relatively small pro- 
portion of the total heat gained by conduction, the resulting error in 
the conduction heat-transfer rate is approximately 1.0 percent. Sub- 
traction of the heat gain through the cork from the total heat gain 
gave values for the heat gained by conduction through the wood, from 
which a conduction correction factor was determined. This correction 
factor is a function of Reynolds number, as the heat conducted throtigh 
the wood was first transferred to the wood by convection from the air 
stream. The values of the heat gained by conduction were obtained by 
using the conduction correction factors and these values were subtracted 
from the total heat gain measured with the surface bare or iced to de- 
termine the convective heat transfer to the plate surface. 

Average skin-friction coefficients for the model both with and 
without ice on the svirface were determined from pressxrre-rake data for 
the subsonic Mach number and pressure -altitude conditions. 
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Mass -transfer , heat -transfer, and skin-friction data are generally- 
presented as a function of Reynolds nimiber. In the calculation of 
Reynolds number for this investigation, the characteristic length used 
was the length from the leading edge of the model to the point where 
the various measiirements were taken. In the case of the sublimation 
experiments, the point at which ice thickness was observed with the 
telescope was used and was approximately 1 inch downstream of the cen- 
ter of the copper cold plate. The heat-transfer Reynolds numbers were 
based on the length from the leading edge to the center of the cold 
plate as the heat-transfer determination was the average for the entire 
cold plate. The length from the leading edge of the wood plate to the 
tip of the pressiire rake (approximately l/2 inch downstream of the 
trailing edge of the cold plate) was used for the characteristic length 
for the friction-coefficient data. Density, viscosity, and Reynolds 
numbers were based on the arithmetic mean of the plate or ice surface 
temperature and the free-stream static tenperature . 


RESULTS AND DISCUSSION 

During the process of sublimation of ice in 'an air- stream, the heat 
required to cause the change of state from ice to vapor is transferred 
to the ice from the air stream. Thus the processes of heat transfer 
and mass transfer are interrelated during sublimation. Experimental 
values of the rate of siiblimation from an ice siirface are used in this 
investigation in conduction with experimental values. of the rate of 
heat transfer from an iced siirface under conditions for which there was 
no transfer of mass to determine values of the Stanton number of sub- 
limation and the ratio of the Stanton nun±>er of sublimation to the 
Stanton number of heat transfer. 

Visual observations of the ice surface indicated differences in 
the ice roughness as a function of (l) the Mach number at which the ice 
was formed, (2) the pressure altitude and Mach number at which subli- 
mation took place, and (3) the time interval from start of sublimation. 
The ice formed at the higher Mach n\mibers was observed to be smoother 
than that formed at the lower Mach number conditions. During the sub- 
limation process, local areas of the ice sTxrface sublimated more rapidly 
than other areas, resulting in a "dished" appearance of the surface; 
however, smoothing out of the individual small projections did occur. 

The dishing effect was particularly severe at the higher-static-pressure 
(low altitude) and low Mach number conditions. At the high-altitude 
condition and particularly at high Mach numbers , the ice became smoother 
during sublimation than it was at the start of sublimation. Since con- 
stant ice thickness was maintained during the heat-transfer study, the 
ice was not subjected to the effects of sublimation on ice-surface 
roughness. 
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Because of the differences in the ice roughness which existed, a 
strict comparison of the sublimation and heat-transfer res\xlts is possi- 
ble only for surface conditions of equivalent roughness. 


Skin-Friction Coefficients 

In order to determine points of equivalent roiighness, a study of 
the sklji friction at siibsonic Mach numbers was made for the sublimation 
and heat-transfer conditions used in this investigation. The study 
also included the determination of skin friction for the bare plate for 
the pvirpose of evaluating the effect of the ice formations on heat 
transfer and skin friction. The results of the skin-friction studies 
are presented in figvire 4 for the bare plate and for the iced plate 
during the subli m ation process. 

The skin-friction data are an average, as equation (lO) has the 
length from the leading edge of the model to the pressure rake as a 
divisor. Thus the skin-friction resvilts are average values existing 
over the model at the instant the data were recorded. The contribution 
to the skin-friction coefficient of that part of the model ahead of the 
cold plate is constant for a given set of Mach number and altitude con- 
ditions. Therefore a direct con 5 )arison may be made of the skin-friction 
coefficients obtained with ice on the surface and with the bare plate 
to determine the effect of the ice formation on skin friction. 

The average skin-friction coefficients obtained by the equation of 
Schoenherr (ref. 4) for a smooth plate and by Baines (ref. 4) for a 
ro-ugh surface are presented for comparison. The curve of Baines pre- 
sented in fig-ure 4 is for x/k = 1500, where k is the roughness pro- 
jection height and x is the surface distance measvired from the lead- 
ing edge. The slope and general shape of the data curves are comparable 
with the CTjrves of Shoenherr and Baines. The effect of Mach nturiber, 
sublimation time, and pressure altitude (hi^er Reynolds mmibers cor- 
responding to low altitudes) on the ro\3ghness of the ice formation and 
consequently on the skin-friction coefficient is indicated. The change 
in average skin-friction coefficient during sublimation. is shown in fig- 
Txre 4 for the 0.4 and 0.8 Mach n\mriber conditions. The tipper dashed 
ctirve for each Mach ntimber shows the skin-friction coefficient shortly 
after the subl im ation process started, and the lower solid curves show 
the skin-friction coefficient after sublimation had progressed for sev- 
eral minutes. The valvies of average skin-friction coefficient show 
negligible change with time at the higher Reynolds numbers for each Mach 
number and appreciable reduction at the low Reynolds numbers, demon- 
strating the smoothing effect of sublimation at high altitude on ice 
roughness. 

In figure 5 is shown the variation of average skin friction with 
Reynolds number obtained with ice of constant thickness (sublimation 
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and condensation occurring at the same rate) on the plate for two Mach 
number conditions. The process of regulating the humidity of the air 
stream to cause the frost point to correspond to the surface tempera- 
ture required a considerable period of time, particularly as the frost 
point approached the surface temperature value. Thus the ice formation 
was subjected to relatively long periods of very slow sublimation or 
condensation, which had the effect of reducing the original ice rough- 
ness. This effect is reflected in the generally lower values of average 
skin-friction coefficient in figure 5 than are shown in figure 4. 

As has been previously noted in the ANALYSIS section, the inclusion 
of the pressure-gradient term in the calculations of average skin- 
friction coefficient would have the effect of increasing the average 
skin-friction coefficient, especially at the higher Mach numbers. Thus 
the curves of skin-friction coefficient in figures 4 and 5 would be 
raised for both the bare- and iced-plate conditions. The bare-plate data 
would be in closer agreement with the Schoenherr line, and the sublima- 
tion skin-friction data (fig. 4) would be brought closer together with 
the 0.6 and 0.8 Mach n\miber lines being raised more than the 0.4 Ifech 
number line. In like manner, the curves of skin— friction coefficient 
in figure 5 would be increased, especially at the higher Reynolds n\mi- 
bers, and the constant Mach number lines would be brought closer to- 
gether. Thus the variation of skin-friction coefficient of the iced 
s“urface with Mach number would be reduced over that shown in the figure. 

Previous measurements (unpublished data) on a dimensionally similar 
flat plate show that the contribution to the local skin-friction coeffi- 
cient by the pres sure -gradient term is of the order of 10 percent for 
the 0.4 Mach number condition and as much as 30 percent at a Mach n-umber 
of 0.8. Application of this type of correction to the data of figures 
4 and 5 cannot be justified because of the lack of boundary- layer data 
with ice on the plate. 


Stanton Numbers of Sublimation and Heat Transfer 

The sublimation data obtained in this investigation were first cal- 
culated by using a value of unity for the ratio Stg/Sth. Values of 
Stg were obtained from equations (2) and (3). The heat-transfer data 
were calculated by using eq-uation (4) to determine Sth. The ratio of 
the Stanton n\msibers thus determined from the experimental data had a 
value of 0.94. Recalculation of the data was then made by using an 
assumed value of 0.90 for the ratio Stg/ Sth wi’th a consequent reduc- 
tion in the values of Stg. The experimental Stanton number ratio was 
then found to be approximately 0.90, which agrees with the assumed value. 
This determination of Stanton nuiriber ratio was made over the Mach num- 
ber range of 0.4 to 0.8 and a pressure-altitude range of 5000 to 30,000 
feet because both these factors influence the ice roiighness. 
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The sublimation data obtained with the value of 0.90 for Stg/Stj^ 

6ire presented in figure 6 in the form of the Stanton number of sublima- 
tion as a function of Reynolds number. The solid faired line was ob- 
tained by use of the method of least squares. The heat-transfer Stanton 
nmber, as a function of Reynolds nxamber, calculated for the case of a 
smooth flat plate from an en 5 >irical equation given in reference 10 is in- 
cluded for comparison. The data show that the Stanton number of subli- 
mation Sts decreases with increasing Reynolds number and has a slope 
very nearly equal to that of the en 5 )irical relation for the heat transfer 
from a smooth plate. No consistent tendency to deviate from the sub- 
sonic resiilts is exhibited by the limited data obtained at a Mach number 
of 1.3. 

Roughness trends as a function of Mach number in the sublimation 
data, which might be expected from the skin-friction results, are not 
apparent in figure 6, probably because of the time -period average factor 
in the determinations of the Stanton ntmibers (roughness changes during 
the period masking the original roughness differences). Another possible 
explanation is that the effect of any ice roughness is of such magnitude 
that changes in the roughness have little effect on sublimation rate. 

In figure 7 is shown the variation of the Stanton number of heat 
transfer Sth as a function of Reynolds nmber for the model both with 
and without ice on the surface of the cold plate. As has been previously 
mentioned in the ANALYSIS section, the values of Sth do not depend on 
assumed values of the Stanton number ratio or the solution of equation 
(2). The eii 5 )irical relation for the Stanton number of heat transfer as 
a function of Reynolds number for a smooth flat plate (ref. lO) is again 
included for con^iarison. The bare-plate data show very good agreement 
with the empirical relation of reference 10. The faired line for the 
iced-plate data was obtained by the method of least squares. Here, as 
was the case with Stg, no Mach number trend was found. Comparison of 

figures 6 and 7 both with and without ice on the surface of the plate 
shows that, for ice roughness of the order of that present during this 
investigation, the sublimation and heat-transfer Stanton numbers are 40 
to 50 percent higher than those for a smooth plate. The increase in 
heat-transfer Stanton numbers due to. ice on the surface was comparable 
with the increase in average skin-friction coefficients caused by ice 
on the surface. In order to determine points of equivalent roughness 
for which true values of the ratio Stg/Sth may be obtained, a plot of 
the skin-friction against Reynolds niomber data presented in figures 4 
and 5 was made as shown in figure 8. These data correspond to the con- 
dition for which the sublimation and heat-transfer data of figures 6 and 
7 were obtained. At the point of intersection of each of the pairs of 
constant Mach number lines, the skin-friction coefficient for the two 
cases is equal and consequently the point of intersection corresponds 
to a point of equivalent roughness. Comparing the value of Stg to 
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that of Stji at the Reynolds numbers corresponding to these two points 
of equal roughness results in values of the ratio Stg/sth of 0.895 and 
0.902 for the 0.6 and 0.8 Mach numiber data, respectively. The average 
value of the experimentally determined ratio Stg/sth is thus approxi- 
mately 0.90, which is in agreement with the value assruned for the cal- 
ctilations. Since the values of Stg and Stj^ of figures 6 and 7 are 

not influenced in any consistent manner by the differences in ice rough- 
ness which were shown to exist by the skin-friction results, a direct 
comparison of the Stanton numbers may be justified over the entire 
Reynolds number range. This ratio of the Stanton numbers is presented 
in figure 9 and has an average value over the entire Reynolds number 
range of approximately 0.90, as was found for the two points of equiva- 
lent roughness. 


Application to Aircraft Icing 

Sublimation of an ice formation from an aircraft surface occ\irs 
during flight in clear air even though the air may be saturated. The 
relatively high surface temperatures of the ice (brovight about by the 
frictional heating associated with high-speed fli^t) raise the vapor 
pressure at the surface of the ice such that it is always greater than 
the vapor pressure of the ambient air. Cvirves showing the sublimation 
rate as a function of altitude and of air velocity over the iced sur- 
faces for two ambient-air temperatures at each of three altitudes are 
shown in figure 10. The value of 0.90 for the Stanton number ratio was 
used in the conputations for these curves. The curves for flight in 
satTjrated air are shown in figure 10(a) and for flight in dry air in 
f igure 10(b ) . A sample problem illustrating the method of calculation 
and the use of the curves of figure 10 is presented in the appendix. 

A comparison of the curves presented in figure 10 shows that the 
rate of sviblimation is somewhat greater for flight in dry air than for 
flight in saturated air, especially for the lower -altitude conditions. 

For example, for a velocity of 550 miles per hovir at a pressure altitude 
of 20,000 feet and an air tenperature of -15° F, 'the rate of sublimation 
in dry air is approximately 15 percent greater than that in saturated 
air; whereas for the 30,000-foot condition at -40° F the increase is 
only approximately 8.5 percent. The decreased difference for the 30,000- 
foot condition results from the very low vapor pressures which accompany 
the low temperatvires at high altitudes. Figiare 10 shows that sublima- 
tion at high altitudes is too slow to be of much value in de-icing of 
aircraft of the transport type which may accumulate ice formations of 
appreciable thickness. 

The dashed curves are presented to show the effect of Reynolds num- 
ber on the rate of sublimation. These curves show that, for the same 
air velocity over the iced surface, sublimation will take place more 
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rapidly nearer the leading edge than at a point some distance back on the 
body. The cvocves are applicable^ however, only for tixrbulent flow over 
the iced surface. 

The curves presented in figure 10 do not necessarily represent the 
maximum sublimation rates which mi^t be experienced in flight. An ex- 
tremely rough ice surface might be expected to yield rates of sublima- 
tion greater than those presented in the figure. 


SUMMARY OF RESULTS 

The experimental investigation to determine the mass transfer by 
sublimation, heat transfer, and skin friction for an iced flat surface 
yielded the following resialts: 

1. The Stanton numbers of sublimation and heat transfer were 40 to 
50 percent higher for an iced s\orface than the Stanton number of heat 
transfer for a bare plate. 

2. For iced surfaces of equivalent roughness, the ratio of the 
Stanton numbers of sublimation and heat transfer was 0.90. 

3. Average skin-friction coefficients obtained with iced siorfaces 
and with a bare plate showed increases due to ice roughness comparable 
with the increase in heat -transfer Stanton number due to ice on the 
surface . 

4. The data obtained at a Mach number of 1.3 showed no appreciable 
deviation in sublimation Stanton number from that obtained at subsonic 
speeds. 


5. Sublimation as a means of removing an ice formation of appre- 
ciable thickness is usually too slow to be of much value in aircraft 
de-icing at high altitudes. 


Lewis Flight Propulsion Laboratory 

National Advisory Committee for Aeronautics 
Cleveland, Ohio, November 18, 1953 
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APPENDIX - SAMPLE CALCULATIONS 

The effectiveness of sublimation as a means of removing an ice 
formation is illustrated by the following numerical example: 

Given conditions: 


Altitude, ft 30,000 

Pressure, in. Hg abs 8.89 

Ambient air temperature, °F -40 

Air velocity over iced surface, ft/sec 805 

Mach number 0.8 

Ice deposit, in 0.2 

Ice specific gravity 0.8 

Ice density, Ib/cu ft 49.9 

Distance from leading edge, ft 1 

Htunidity, at ambient air temperature saturated 

Latent heat of sublimation, Btu/lb 1220 

Specific heat of air, Btu/(lb)(°F) 0.24 


Ratio of molecular weight of water to molecular weight of air. . .0.622 

Problem: To determine the time required to remove the ice by 

sublimation. 

Equation (2) is used to determine the surface temperature and sur- 
face vapor pressuire of the ice, and equation (3) is used to determine 
the weight rate of sublimation. The value of 0.90 will be used for the 
ratio Sts/Sth, and the value of Sts is taken from figure 6. 

Equation (2) may be. restated as 


t -+ ( ^s ^1 \^s 

S ■ St^ m^^Pl - eg ■ Pi - e^^Cp 

The adiabatic wall teniperatiire is obtained from the equation 

tad = ^o(l + 0.176m2) 


( 2 ) 


where tg^^i and tQ are in °R. Thus tg^^ = 467.3° R or 7.3° F. The 

saturated vapor pressure at the ambient temperat\ire is 0.0039 inch 
mercury from reference 11. Substitution of the known values into the 
equation for determining surface temperature results in the following 
equation : 
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ts 


7.3 - 0.90 X 0.622 X 


1220 / 

0.24 \8.89 - eg 


0.0039 

8.89 - 0.0039 


) 


By assviming values of ts and corresponding values of es from ref- 
erence 11 and solving by trial-and-error procedure, the value of ts 
foxind to be -2.3° F with a corresponding value of es = 0.0338 inch 


mercury. The q\iantity 


vPl - 



has a value of 0.003378, 


IS 


The ten 5 )erature which is used in determining the air viscosity, air 
density, and Reynolds number is the arithmetric mean of the ambient air 
temperature and the surface teit 5 >erature ; or (tQ + tg)/2 = -21.2° F. 

The air density is fovind to be 0.0269 pound per cubic foot and the 
Reynolds number at a point 1 foot from the leading edge is 2.066X10°. 
From figure 6 the value of Stg is fomd to be 2.46X10"^. 


Equation (3) solved for the weight rate of sublimation becomes 


W = StgpU 


^ / es 
®a \Pl - °s 



Substituting values gives 

W = 2.46X10"^ X 0.0269 X 805 X 0.622 X 0.003378 
= 0.000112 lb/(sec)(sq ft) 

= 0.403 lb/(hr)(sq ft) 

The weight of ice to be removed from each sqiiare foot of s-urface is 
Ice weight = 0.20/l2 ft X 1 sq ft x 49.9 lb ice/cu ft 
= 0.832 lb 


The time for removing the ice hy suhlimation is thus 


0.832 Ib/sq ft 
0.403 lb/(hr)(sq ft) 


2.06 hr 


The solution of a problem with conditions within the range of fig- 
\xce 10 may be readily obtained as demonstrated by using the conditions 
of the example. For a velocity over the iced surface of 805 feet per 
second (549 raph) , the rate of sublimation given by figure 10(a) is again 
0.403 pound per hour per sq\aare foot. 
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For the case where sublimation takes place in completely dry air, 
figinre 10(b) yields a sublimation rate of 0.432 pound per hour per 
sqviare foot; and the time for removal is 1.93 hoiirs. 

If atmospheric conditions were such that the ice formation in the 
preceding exan 5 )le could be sublimated at an altitude of 20,000 feet and 
at an air temperature of -5° F, the sublimation rate as found from fig- 
ure 10(a) for saturated air is 1.103 pounds per hour per square foot 
and the removal time is approximately 45 minutes. 
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Figure 2. - Schematic diagram of apparatus for suhllmatlon, heat-transfer, and skln-frlctlon studies of Ice 
on a flat plate. 
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Figure 3, - Photograph of Ic© foniatlon on surface of cold plate shoving typical presuhlimation condition 
of iced surface. 


Time after start 
of sublimation. 
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Figure 5. - Experimental values of average skin- friction coefficient for iced surface at 
constant ice thickness with heat transfer from air to ice. 


Stanton number of heat Stanton number of sublimation 

transfer, St^ Stg 
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Figiare 6. - Variation of sublimation Stanton number with Reynolds 
number for four values of Mach number. Assiomed Stg/Stj^, 0.90. 



Figure 7. - Variation of heat-transfer Stanton number with Reynolds 
number for bare plate and iced surface. 





Velocity over Iced surface, mph 

(a) Saturated air at ambient temperature. 

Figure 10. - Sublimation rate as function of velocity over 
Iced surface for two values of temperature at each of three 
altitudes . 








